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We use the Renormalization Group method to study the magnetic field influence on the Bose- 
Einstein condensation of interacting dilute magnons in three dimensional spin systems. We first 
considered a model with SU(2) symmetry (universality class z — 1) and we obtain for the critical 
magnetic field a power law dependence on the critical temperature, [H C (T) — H c (0)] ~ T 2 . In the 
case of U(l) symmetry (universality class z = 2) the dependence is different, and the magnetic 
critical field depends linearly on the critical temperature, [H C (T) — H C (G)] ~ T. By considering a 
more relevant model, which includes also the system's anisotropy, we obtain for the same symmetry 
class a T 3,/2 dependence of the magnetic critical field on the critical temperature. We discuss these 
theoretical predictions of the renormalization group in connection with experimental results reported 
in the literature. 



I. INTRODUCTION 

Bose-Einstein condensation (BEC) remains one of the 
most important macroscopic effects predicted long time 
ago by quantum mechanics for an ideal Bose gas. Lately, 
the interest in BEC was renewed by the practical real- 
ization of a condensate phase in weak interacting Bose 
systems realized in ultra-cooled diluted alkali atomic 
gasesi. Further, a BEC of fermionic pairs was achieved 
in trapped fermionic systems under attractive fermion- 
fermion interaction^. The limitation of these experi- 
ments is due to a reduced number of particles which 
undergo BEC. Recently, BEC was associated with the 
magnetic transition observed in different quantum spin 
systems. A particular class of such materials is formed 
by the XCuCi3 dimer compounds (with X being Tl, K, or 
NH,i)2ii. In the ground state these materials are spin sin- 
glets which are magnetically inactive. If a high-magnetic 
field is applied, they undergo a transition into a magneti- 
cally ordered state, a transition which can be understood 
as a condensation of excitations which behaves as bosonic 
quasiparticles. 

One particular example is TICUCI3. The quantum 
magnetism in this compound is attributed to the spin 1/2 
Cu 2+ ions positioned in the double chains of CU2CI6. The 
ground state of TICUCI3 was found to be a singlet with an 
excitation gap A ~ 7.5 K£. The gap has been associated 
with the weak anisotropic antiferromagnetic (AF) intra- 
dimer coupling in the double chain. Magnetic suscep- 
tibility experiments at different temperatures for differ- 
ent directions of the applied magnetic field exhibit broad 
maxima at T = 38 K and a decreasing behavior toward 
zero susceptibility with decreasing temperature^. At very 
low temperatures, the magnetic susceptibility for differ- 
ent values of the applied magnetic field behaves quite 
differently. At fields of the order of 1 T, the magnetic 
susceptibility decreases exponentially to zero, a proof 
for the existence of a ground state gap. On the other 
hand, for a magnetic field of 7 T an anomaly in the 



magnetic susceptibility was reported around T = 4 K. 
This is an indication that the ground state must be gap- 
less at high magnetic fields. The system has a three di- 
mensional character and in an external magnetic field 
H, the singlet-triplet gap A is reduced to A — fisgH, 
and vanishes at a critical field H c = A/[1b9- Inelastic 
neutron scattering measurements proved the existence 
of the elementary magnonic excitations with a strong 
dispersion in all three directions&t The observed gap, 
A = 0.7 meV, which is much smaller than the intra- 
dimer interaction J — 5 meV, and the small critical field 
H c = 5.6 T, make this material a perfect candidate to 
study the magnetic field induced phase transition. An- 
other argument in favor of the possible phase transition 
was the evidence for a Goldstone mode& observed also 
by inelastic neutron scattering experiments^. Neutron 
diffraction experiments at fields H > H c showed that 
the field induced AF order in the plane normal to the 
applied field appears at the same time with the uniform 
moments^. Sound attenuation experiments were associ- 
ated to a "relativistic"-like form of the excitations energy 
spectrum, E(k) = ^/A 2 + k 2 /[2m e ff] (m e ff — 2/ J 2 ), a 
relation which was used to explain the connection be- 
tween the transition temperature and the concentration 
of magnonsii. The main conclusion of all these experi- 
mental approaches is that TICUCI3 undergoes a magnetic 
ordering in high magnetic fields. The temperature depen- 
dence of the critical magnetic field can be summarized by 
[H C (T) — H C (Q)] oc T^, where from the experimental data 
fit <j)~ 1.7-2.2. 

Various theoretical models were used to investigate the 
BEC of magnons in quantum spin liquids. In the mean 
field approximation most of the experimental features of 
the phase transition cannot be reproduced, in particular, 
the critical temperature dependence on the magnetic field 
being almost fla1pi£. The possibility of a magnetic field in- 
duced phase transition in quantum spin liquids was first 
discussed in terms of BEC by Giamarchi and Tsvclik 13 . 
This possibility was investigated using the Hartree-Fock 



2 



approximation introduced by Popoxii*. The initial sys- 
tem is mapped into a dilute Bose liquid with a magnetic 
field dependent chemical potential \x = guB^H — H c ), 
and the total number of magnons N is associated with 
the total magnetization M, M — gfi B N—. The theory 
predicts [H C {T) - H c (0)} ~ T 3 / 2 , a result which is not in 
complete agreement with experimental data. The model 
was further on investigated by Misguich and Oshikawa 16 
by including a more realistic dispersion relation^ and the 
shape of the critical temperature curves as function of the 
critical magnetic field was well reproduced with a single 
adjustable parameter. 

The possible phase transition was numerically inves- 
tigated by Monte Carlo simulationiLiLiS. Nohadani 
et al*L found that the value of the critical tempera- 
ture exponent <j> is in general greater that the Hartree- 
Fock value <fi — 1-5, but is close to it as the mag- 
netic field approaches its critical zero temperature value. 
Kawashima 18 argued that the Hartree-Fock value <j> = 
1.5 may be incorrect at finite temperature, but is the 
correct value in the T = K limit, namely for the sit- 
uation when the phase transition is actually a quantum 
phase transition (QPT). The BEC of magnons as a QPT 
at T = K was further investigated by Monte Carlo 
simulation by Nohadani et al. 19 . The model considered 
a control parameter g — J'/ J, where J is the intra- 
dimer exchange interaction, and J' is the inter-dimer 
coupling. The experimental data gives for these cou- 
plings the values J = 60 - 70 K and J' = 40 - 53 K. 
However, those values can be controled by external pres- 
sure. The resulting T = K phase diagram consists of 
three different regions, regulated by the g ratio. At low 
magnetic fields and small control parameter g the system 
is a dimer spin liquid, i.e., is in a magnetically disordered 
state. This phase is characterized by a SU (2) symmetry 
with a dynamical critical exponent z = 1. At intermedi- 
ate magnetic fields H < H s (H s is the saturation field) 
and sufficiently large value of g, the ground state is par- 
tially spin polarized for H > H c (H c is the critical field) 
and has a long-range antiferromagnetic order transverse 
to the applied magnetic field direction. In this case the 
phase symmetry is only 17(1) with a corresponding dy- 
namical critical exponent z = 2, a value associated to 
the quadratic dispersion of the bosonic excitations. At 
large fields H > H s (H s ~ J + 5 J'), all spins are fully 
polarized. 

Here we consider a Renormalization Group (RG) ap- 
proach of the QPT associated to the BEC of magnons in 
spin dimer systems. It is well known that the presence 
of a QPT will influence the system's properties even at 
finite temperatures. The exact temperature range in the 
vicinity of the T = K phase transition is unknown, 
however, as in any phase transition a critical region will 
be present. We will consider separately the phases identi- 
fied by Monte Carlo simulation, and show that different 
values for the critical temperature exponent <j> are ob- 
tained in the RG analysis. We will show that in the 
phase with a SU(2) symmetry cj> = 2, whereas for the 



phase with J7(l) symmetry <f> — 1. If one considers a 
more realistic model with a U{\) symmetry and strong 
anisotropy, (p = 3/2. This behavior suggests the exis- 
tence of a narrow crossover temperature interval close to 
T = K (T C 1 K) where the critical temperature co- 
efficient complies with the Hartree-Fock value = 3/2. 
At higher temperature, still in the critical region, the 
symmetry of the system changes and the value of the 
temperature critical exponent is higher <fi ~ 2. We will 
discuss those results in connection with the numerical 
calculations and experimental data. 



II. RENORMALIZATION GROUP ANALYSIS 

Following Herta 2 ^, and MilliaSi, we describe the critical 
behavior of a magnon dilute gas in terms of an effective 
Ginzburg-Landau- Wilson theory of the order parame- 
ter field $(k, t) which represents the fluctuations of the 
staggered magnetization of the system. 

In the general theory of QPT's the total action will 
consist on two terms, S[$] = S 2 [§} + 5" 4 [$]. The 
quadratic part, S , 2[$], takes the form 

= \ E ^(k)[r (H) + k° + (iu n r}<t>(k) , (1) 

fe 

where k = (k, uj n ) with u> n — 2nnT (fcs = 1), <&(k) is 
a bosonic field describing the magnetization fluctuations, 
and ro (H) measures the distance from the quantum crit- 
ical point. In the presence of an external magnetic field 
the control parameter r (H) will acquire a field depen- 
dence. The form of the action was extensively discussed 
by Fisher and Rosch 2 ^ both on a phenomenological ba- 
sis and starting from a Hubbard-type model of electrons 
moving in the presence of a magnetic field. In the follow- 
ing we will consider to have a Zeeman type H dependence 
of the control parameter, tq(H) ~ [H — H c ]. The dimen- 
sionality of the system is d = 3, and a and m will take 
different values according to the studied phase, i.e., a = 2 
and m = 2 meaning that the dynamical critical exponent 
z = 1 for the SU (2) symmetry, and a = 2 and m = 1 
with a dynamical critical exponent z = 2 for the U{\) 
symmetry. The interacting contribution to the action, 
54 [$] , will be of standard form 

St[«] = 

^ E *(fci)*(fe)$ t (fc 3 )$ t (fc4)<5(fc 1 + fc 2 + fc 3 + fc 4 ) , 
16 

k i , . . . , fe'4 



and describes interactions between the considered fluc- 
tuations. The total action will remain invariant under 
a standard transformation of momenta, frequency, and 
fields, i.e., k' = kb, uj' n = u) n b z , where z is the dynamical 
critical exponent, and $' = $6~( d + z + 2 )/ 2 M 

The scaling equations for the parameters T, r, and u 
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have the general form 24,26 : 



= zT(l) , 

2r(l)+K 3 F 1 [r{l),T(l)] , 



(3) 
(4) 



dT{l) 

dl 
dr(l) 
dl 

^ = [4- (d + z)]u(l)- b -K 3 F 2 [r{l),T{l)]u 2 {l) ,(5) 

where Fi[r(l),T(l)] and F 2 [r(l),T(l)} are complicated 
functions of parameters T and r. Their values in the 
low temperature domain and at the critical point (7' = 0) 
are given in Ref£ 4 as Fi \T(l)] = coth [l/(2T(Z))]/2 and 
F2[T(l)] — 1/4. Also, K 3 is a constant whose value is 
l/(27r 2 ). Equations ©-I© will be solved to analyse the 
system behavior in the critical region for two situations. 
First we will analyse the SU(2) symmetry situation when 
the dynamical critical exponent z = 1, and then the U (1) 
symmetry case when z = 2. 



where t r (0) = ro — tq c with tq c — —K 3 uq/4, and ((x) is 
the Riemann zeta function. The renormalization proce- 
dure will be stopped at / = I* ^> 1, given by t r (l*) = 1. 
To find I* we rewrite t r (l) — e 2l t r (T) and after simple 
calculations one finds 



I* 



ml 

T 



(12) 



We can define a critical line whose equation is given by 
t r (T) = which permits us to calculate the critical field 
H C (T) as: 



H C (T) = H c {0) - K 3 u C,{2)T 2 



(13) 



Accordingly, in the SU(2) symmetry phase we can iden- 
tify the exponent of the critical temperature dependence 
on the magnetic field to be = 2, a value which is close 
to the one reported in the literature. 



A. SU(2) symmetry 

One possible phase identified by Monte Carlo simula- 
tions was characterized by the limit of low magnetic fields 
and small coupling ratios J' /J, being a magnetically dis- 
order phase, i.e., a dimer spin liquid. In this case the 
dynamical critical coefficient value is z = 1. Two of the 
renormalization group equations, namely © and JSJ, can 
be solved exactly with the following solutions 



T(l) = T e l 



and 



u(l) 



C (l + l ) 



(6) 



(7) 



where Co = 5K 3 /8 and Iq = l/[Co«o]- The remain- 
ing equation the one for the control parameter r, 
will be solved considering a solution of the form r(l) ~ 
ro exp(2l)h(l), leading to 



r(l) 



where 



K 3 u(l) 



1 dx exp[—2x\u{x) 
exp [l/T{x)\ - 1 



(8) 



(9) 



To study the influence of the criticality on the thermody- 
namics of the SU (2) phase we introduce the scaling field 
t r (l) defined as: 



*r(0 = r(0 



K, 



u(l) 



(10) 



a field which will be used to stop the renormalization 
process. Based on Eqs. Q and JSJ we get 



,,2/ 



"*r(0) + K 3 a2)u T 2 



(11) 



B. U(l)symmetry 

Another possible phase identified by Monte Carlo sim- 
ulations of the quantum phase transition in spin dimer 
liquids has U(l) symmetry and is present for intermedi- 
ate applied magnetic fields and a sufficiently large value 
of the ratio J'/ J. In this case the dynamical critical 
exponent value is z — 2, and the renormalization group 
equations ©-© are similar to those of the weakly inter- 
acting Bose gas^S&. Equation J2J has a trivial solution 
of the form 



T(l) = Te 



21 



(14) 



We consider now equation (JjJJ which gives us the renor- 
malization of the interaction parameter. If we introduce 
e = 2 — d < (d=3) this equation can be rewritten in 
the form 



du(l) 



eu(l) 



K \ 2 (l) 



dl 4 
and admits the following solution 

M exp[-|e|;] 



u(l) 



l + X 3 (exp[-|e|;]-l)/4 



(15) 



(16) 



One can proceed now to the solution for the renormalized 
control parameter r. In the U(l) symmetry case this 
equation will admit the following solution 



r(i) 



K 3 u(l) 



K 3 u , K 3 Tu{l) 



hi 



1 



1 - cxp[-l/T(0] 



(17) 



As in the previous situation we introduce the scaling field 
t r (l) — r(l) + K 3 u{l) I '4 and we will stop the renormaliza- 
tion procedure at t r (l*) = 1. With this definition, based 
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on equation (|17(l . the general form of the scaling field 
become 



t r (l) = e 21 { t r (0) 



K 3 u{l) 



Tin 



1 



1 - exp[-l/T(0] 



(18) 

where t r (0) — r$ — r$ c with vq c = —K 3 uo/A. We stop 
the scaling procedure at I = I* 3> 1, where once again 
I* is the solution of t r (l*) = 1. In this case the equation 
is more complicated, however, in the low temperature 
regime we can approximate its solution. If we rewrite 
t r (l) = e 2l t r {T) with 



t r {T) ~ r - r , 
we finally find 

r - 



1 

——w In — 

2 uo 



2 T 



(20) 



The critical line, identified from t r (T) = will give us the 
critical magnetic field H C (T) as function of temperature 



H C (T) = H c (0) 



uo , 1 

47T 2 



u 



T. 



(21) 



In the case of U(l) symmetry the relation between the 
critical field and temperature is different from the case of 
SU(2) symmetry and we identified in this situation the 
critical temperature exponent to be <j) = 1. 



dr{l) 
dl 



= Ml) 



-^{[v(l) + 2u(l)]F 1 [T(l)} + v(l)T(l)},(2A) 



du{l) 

dl 
dv(l) 

dl 



:eu(0-^-V(0 
ev(l) . 



(25) 
(26) 



The last two equations, the ones corresponding to the 
renormalized coupling constants u(l) and v(l), have been 
written in the lowest order, an approximation which is 
assumed valid in the low temperature domain. Equation 
(19) (HSJl admits the trivial solution T(Z) = Te 21 . The first 
coupling constant, u(l), has the same form no matter 
if one considers or not the anisotropy, given by equa- 
tion i|16|) . Equation l|2t)[l admits the trivial solution 



v(l) 



v e 



\el 



The solution for the remaining equa- 



tion which is giving us the renormalized phase transition 
control parameter, r, can be written as 



r(l) = e 2 



^L^dl'v(l')T(l') 



' 2 . hl 



dl'[v(l')+2u(l')}F 1 [T(l')}},(27) 



where ro c = K^vq + 2uq]/8. We can introduce again the 
scaling field t r (l) = r(l) + K 3 [v(l) + 2u(l)}/8 which based 
on the solution for r can be written as 



III. INFLUENCE OF ANISOTROPY 

According to Monte Carlo calculations at sufficiently 
large fields the system has a ground state which is par- 
tially spin polarized and has antifcrromagnetic long- 
range order transverse to the magnetic field direction. 
The fact that the magnetic properties of the system de- 
pend on direction suggests that anisotropy may play an 
important role when the system is in this partially po- 
larized phase. To take in account such effects, we con- 
sider that the interaction term in the total action corre- 
sponding to the 17(1) symmetry phase has the following 
formM 



^ = JQ E $(^(fe)^3)^fc4)<5k 1+k2 ;k 3+ k 4 

k± . . . .,&4 

(22) 

where the coupling constants uq and vq describe interac- 
tions between the magnetization fluctuations in different 
directions of the system. Using the same procedure as in 
the previous cases the renormalization group equations 
can be obtained as 



dT(l) 
dl 



2T(l), 



(23) 



U(l) 



e 21 < r 



roc 



K 3 v Q 
2e 



T[e- |e|i - 1] 



K 3r 



+^-T[v(l) + 2u(l)] In 



1 



1 - exp[-l/T(Z)] 



(28) 



where e < for d = 3. To calculate the value of the stop 
scaling parameter I* we consider t r (l*) = 1 and, in the 
lowest order, we obtain I* as: 



I* 



iml. 

2 T 



(29) 



Equation l|28|l for the scaling field can be rewritten in the 
form t r (l) = e 2l t r (T) where t r (T) can be evaluated as: 



tr(T) 



K 3 v T 
2\e\ 



+ ^( Vo + 2uo)T"/ 2 ln 



1 



v + 2u 



(30) 



In the anisotropic U(l) case the temperature dependence 
of the critical magnetic field will be calculated in the limit 
«o > vq from t r (T) = as 



H C (T) ~ H c (0) - 



uo , 1 
— 7T In — 

2tt 2 



2u Q 



T 3/2 



(31) 
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A similar result was predicted using different calculation 
methods and numerical evaluation in Refs^iS. Equa- 
tion predicts an important role for the anisotropy in 
the case of U(l) symmetry, as the critical temperature 
exponent changes from <j> = 1 to = 3/2. The later 
value is close to the lowest observed experimental value, 
i.e. 1.7, and identical to the value reported by Monte 
Carlo studies. 



IV. DISCUSSIONS 

The idea of BEC in solid state systems was associated 
with electronic Cooper pairs in superconductors, excitons 
in semiconductors, and more recently with magnons in 
spin liquid dimer compounds. The occurrence of a con- 
densate phase in these compounds was investigated by 
various theoretical methods, from mean field to renor- 
malization group approximations. In the case of spin liq- 
uid dimer compounds the BEC of magnons is induced by 
magnetic fields, the critical field associated to the BEC 
being temperature dependent and characterized by a crit- 
ical exponent 0, defined as \H C (T) - H c (0)} ~ T*. The 
exact value of the critical exponent <f> can vary from 1.7 
to 2.2 according to experimental data. A mean field anal- 
ysis gives (f> — 1.5, and for some time it was believed that 
such approximation was enough to explain all relevant 
physics of the phase transition. 

Here we applied the renormalization group method to 
study the possible BEC of magnons. We started our 
analysis from the premise that a quantum phase transi- 
tion, despite being characterized by a T — K critical 
temperature, will influence the system properties even at 
finite temperatures. Accordingly, we analyze the phase 
transition in the vicinity of a quantum phase transition, 
in the low temperature limit. The possible quantum 
phase transition in a spin liquid dimer was investigated 
by Monte Carlo simulations considering a model which 
includes also the system's anisotropjsiS. As function of 
anisotropy the phase diagram of a spin liquid dimer will 
consists of three different regions, in each of them the sys- 
tem being in a different symmetry class At low mag- 
netic fields, the system is in a magnetically disordered 
state. As the fields increases, a partially spin polarized 
state will develop and the symmetry of the system will 
change accordingly. At high magnetic fields the system is 
in a long-range order antiferromagnetic phase. To take 
into account all these possibilities, we analyzed the in- 



fluence of magnetic fields and temperature on the BEC 
of magnons in systems with different symmetry. In the 
case of a magnetically disordered spin liquid dimer, when 
the symmetry of the system is in the SU(2) class and 
the dynamical critical exponent is z = 1, the critical 
exponent is <f> = 2. On the other hand, when the par- 
tially ordered phase is induced by the magnetic field, the 
system's symmetry changes to E/(l), and the dynami- 
cal critical exponent becomes z = 2. In this situation, 
we calculated <f> = 1, a value which was never observed 
in real experiments. However, in this partially ordered 
state anisotropy plays an important role, as the phase 
is characterized by a large J' / J ratio. When anisotropy 
is taken into account by considering direction dependent 
interactions in the action, the value of the critical coef- 
ficient (j) changes from 1 to 1.5, a value which is close to 
the lowest reported experimental data. The influence of 
anisotropy was considered by Fischer and Rosh 22 includ- 
ing additional terms in the free action, the final result 
for the temperature critical exponent being <f> = 1.5, the 
same value we obtained. The crossover problem was also 
considered in a sigma model framework to explain some 
phases of the organic insulator (TMTTF^PFg as func- 
tion of temperature, magnetic field, and pressureiSL 2 ^ A 
similar idea was used to explain nuclear magnetic reso- 
nance experimental data in superconducting spin-ladder 
compounds such as Sr2Cai2Cu2404iiSS 

In conclusion, we showed that the renormalization 
group analysis of the possible BEC of magnons in spin 
liquid dimers is a more appropriate investigation method. 
This method can account for the change in the critical 
exponent <fi from 1.5 to 2, according to the phase symme- 
try. The results are in good agreement with the experi- 
mental data in TICUCI3 samples. More recently a simi- 
lar behavior was reported also in CS2CUCI4 samples^. 
Such a change in the critical exponent cannot be ex- 
plained by mean field approximations. On the other 
hand the ground state of TICUCI3 is strongly influenced 
by pressured Accurate experimental data obtained for 
various values of the applied external pressure still sup- 
port previous results for the critical exponent, namely 
<fi = 2. 6. 31 Finally, we mention that recent experimental 
data showed that the occurrence of the BEC of magnons 
can be also induced by hydrostatic pressure 3 ^. The de- 
pendence of the critical field on pressure, and the pressure 
dependence of the transition temperature are suggesting 
a similar theoretical description. 
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